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Abstract 

A canal surface is an envelope of a one parameter family of spheres. In 
this paper we present an efficient algorithm for computing the implicit 
equation of a canal surface generated by a rational family of spheres. By 
using Laguerre and Lie geometries, we relate the equation of the canal sur- 
face to the equation of a dual variety of a certain curve in 5-dimensional 
projective space. We define the /i-basis for arbitrary dimension and give 
a simple algorithm for its computation. This is then applied to the dual 
variety, which allows us to deduce the implicit equations of the the dual 
variety, the canal surface and any offset to the canal surface. 
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1 Introduction 

In surface design, the user often needs to perform rounding or filleting between 
two intersecting surfaces. Mathematically, the surface used in making the round- 
ing is defined as the envelope of a family of spheres which are tangent to both 
surfaces. This envelope of spheres centered at c(t) £ R 3 with radius r(t), where 
c(t) and r(t) are rational functions, is called a canal surface with spine curve 
£ = {(c(t),r(t)) € R 4 |i S R}. If the radius r(t) is constant the surface is called 
a pipe surface. Moreover, if additionally we reduce the dimension (take c(t) in a 
plane and consider circles instead of spheres) we obtain the offset to the curve. 
Canal surfaces are very popular in Geometric Modelling, as they can be used as 
a blending surface between two surfaces. For example, any two circular cones 
with a common inscribed sphere can be blended by a part of a Dupin cyclide 
bounded by two circles as it was shown by |Pratt(1990)| Pratt(1995)| (see Fig- 



ure[T]). Cyclides are envelopes of special quadratic families of spheres. For other 
examples of blending with canal surfaces we refer to |Kazakeviciute(2005)| . 

Here we study the implicit equation of a canal surface C and its implicit de- 
gree. The implicit equation of a canal surface can be obtained after elimination 
of the family variable t from the system of two equations gi(y, t) — g?,{y, t) = 
(here 91,(72 are quadratic in the variables y — (3/1,2/2,2/3,2/4)), i-e. by taking the 
resultant with respect to t. However, this resultant can have extraneous factors. 
In the paper we explain how these factors appear and how we can eliminate 
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Figure 1: A Dupin cyclide used for blending circular cones. 



them. By using Lie and Laguerre geometry, we see that the above system of 
equations is related to a system h\{y,t) = fi2(y,t) = Q(y) = 0, where h%,h2 
are linear in the variables y = (u, yo, y±, y2, V3, Va) and Q(y) is the Lie quadric 
(for the exact definition see formula ©). It turns out that the variety defined 
by the system of equations hi(y, t) = h.2(y, t) = is a dual variety to the curve 
£ G P 5 , where £ is a curve on the Lie quadric determined by the spine curve 
£ (for the explicit definition see formula (flT|) ). For the dual variety V(£) we 
define the /i-basis, which consists of two polynomials Pi{y,t),p2{y,t) which are 
linear in y, and of degree d\,d,2 in t are such that d\ + d,2 is minimal. It turns 
out that the resultant of p\ and P2 with respect to t gives the implicit equation 
of the variety V(£). There is a simple substitution formula (see the algorithm 
at the end of section [5|) to compute the implicit equation of the canal surface 
from the implicit equation of the variety V(£). 

Partial solutions to the problem of finding the implicit equation (and degree) 
for canal surfaces have been given in other papers. For instance, the degree 
of offsets to curves is studied in Se gundo, Sendra(2005)| . In Xu et al.(2006)| , 
there is a degree formula for the implicit equation of a polynomial canal surface. 
Quadratic canal surfaces (parametric and implicit representation) have been 
studied in Krasauskas, Zube(2007)| . 

We close the introduction by noting that the implicit degree of a canal surface 
is important for the parametric degree. Our observation is that if the canal 
surface has the minimal parametrization of bi-degrec (2, d) then its implicit 
degree is close to 2d. On the minimal bi-degree (2,d) parametrizations of the 
canal surface we refer to [Krasauskas(2007)] . 

The paper is organized as follows. In the next section we develop some 
algebraic formalism about modules with two quasi-generators. We define the 
/^-basis for these modules and present an algorithm for its computation. In the 
following section, we recall some needed facts about Lie and Laguerre sphere 
geometry. Then using Lie and Laguerre geometry we describe the canal surface 
explicitly Also, we introduce the T-hypersurface which contains all d-offsets to 
the canal surface. Using the /i-basis algorithm we compute the implicit equations 
of the dual variety V(£), the T-hypersurface and the canal surface C. Next we 
apply the results of the previous section to the dual variety V(£ ) of the curve 
and explain how to compute the implicit degree of the T-hypersurface (without 
computation of the implicit equation). Finally, we give some computational 
examples. 
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2 Modules with two quasi- generators and the //- 
basis. 



Let R[t] be polynomial ring over the field of real numbers, and denote M.[t] d the 
i?-module of <i-dimensional row vectors with entries in R[t]. Let R(t) be the 
field of rational functions in t. For a pair of vectors A — (Ai, A 2 , . . . , A d ),B = 
(B x , B 2 , B d ) eR[t] d the set 

M = (A, B) = {aA + bB e R[t} d \ a,b G R(t), A, B e R[t] d } c R[t] d (1) 

is the R[t]-modulc with two polynomial quasi- generators A, B. Here, we assume 
that A, B arc R[t] -linearly independent, i.e. aA + bB = with a, b G R[t] if and 
only if a = b = 0. 

Remark: Note that the vectors A, B may not be generators of the module M 
over R[t] because a and b in the definition (JTJ) are from the field R(t) of rational 
functions. For example, if A = pD with p G R[£], -D G R[i] d and degp > then 
A, B are not generators of the module M. 

For A = (A 1 ,A 2 ,...,A d ),B = (B 1: B 2 , ■-, B d ) G M[t] d we define the Pliicker 
coordinate vector A A B as follows: 

iAB= ([1,2], [1,3],.., [d- G R[i] d(d - 1)/2 , where = AiBj - AjBi. 
In other words, A A B is the vector of 2-minors of the matrix 

w AB = ( Al M ■ ■ Ad ) 

VVA < B \Bx B 2 ••• B d J 

and we denote by deg(A A B) = max iiJ {deg(A i _B :) - — AjBi)} the degree of the 
Pliicker coordinate vector, i.e. the maximal degree of a 2-minor of Wa,b- 
Let a polynomial vector A G M.[t] d be presented as 

n 

A = V" ait* , a, G M d , i = 0, .., n; a„ ^ 0. 

i=0 

We denote the leading vector a n by LV(j4) and the degree of A by deg^4 = n. 

Note that if LV(A) and LV(B) are linearly independent over R then deg A A 
B = deg A + degB and LV(A A B) = LV{A) A We define 

deg M = min{deg(A A B) A,B G R[t] d such that (A, B) = M} 

to be the degree of the module M with two quasi-generators. 

Definition 1. Two quasi-generators A, B of the module M = (A, B) are called 
a //-basis of the module M if degM = deg A + degB. 

As we always have the inequality deg(A A B) < deg A + deg B, this means in 
particular that the sum deg A + degi? is minimal. A /i-basis always exists, as 
we shall see at the end of the section. Let us explain the geometric motivation 
behind this definition. 
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Remark 2. By abuse of notation, we will continue to denote parameters t, 
however in the geometric definitions that follow, they should be understood 
as parameters (t : s) € P 1 and polynomials in M.[t] should be thought of as 
homogenized with respect to a new variable s. 

We define the following subspace of M. d for the module M = (A, B). 

L(M, t ) = {x e R d | C(t ) • x = for all C 6 M} 

where C(t) = {C x {t), C 2 (t), . . . , C d (t)) £ R[t] d , x = (x u x 2 , . . . , x d ) T and C(t) ■ 
x = x\Ci(t) + x 2 C2(t) + . . . + XdCd{t). We have the inequality dim L(M, t ) > 
d — 2, because the module M has only two quasi-generators. In fact, we have 
dim L(M, to) = d — 2 for all to, as we will see in Proposition [HI 2. Whenever 
two vectors A(to) and B(to) are linearly independent in E d then L(M, to) is the 
intersection of two hyperspaces {n e l^j A(to) -x — 0} and {x 6 M d | B(to) -x = 
0}. 

Using those subspaces, we can associate a hypersurface Sm in the real pro- 
jective space P^ 1 = P(R d ) with the module M 

S M :=\JV(L{M,t)) CP^ 1 . (2) 
t 

Note that this definition and the definition of L(M,to) depend only on the 
module M and not on the choice of quasi-generators. It is useful to compare 
the hypersurface Sm with the hypersurface Sa, b defined as 

S a ,b ■■= \J({A(t) ■ x} n {B(t) ■ x}) c P^ 1 (3) 

t 

where A, B are quasi-generators of M. By definition, this is the variety de- 
fined by Res t(A(t) ■ x, B(t) ■ x) and it is clear that Sm C Sa,b- If the vectors 
A(t ),B(t ) arc linearly dependent, then {{A(t Q ) ■ x} n {B(t Q ) ■ x}) C K d is 
a subspace of codimension one. Note that in this case the implicit equation 
Rest(A(t) ■ x,B(t) ■ x) contains the factor A(to) ■ x. As a matter of fact, this 
happens if and only if WA,B(to) has rank one, which is equivalent to saying that 
to is a zero of the ideal generated by the Pliicker coordinates. 

In fact, we will see in Proposition [6] that this phenomenon does not occur 
for /i-bases, i.e. if A, B is a /u-basis of the module M then Sm = S^ ^ and there 
are no extraneous factors as before. 

Remark 3. We should explain why we use the term /i-basis. The above def- 
inition is a generalization of the usual definition for the ^-basis of a ratio- 
nal ruled surface (as in Cox, Sederberg, Chen (1998)] , [Chen et al.(2001)] or 



Dohm(2006)| ). They coincide in the special case d = 4. M is the analogue of 
the syzygy module (i.e. the module of moving planes following the parametriza- 
tion of the ruled surface) and the subspaces L(M,t), which in this case are 
2-dimensional and hence define projective lines, are exactly the family of lines 
which constitute the ruled surface. Similarly, the case d = 3 corresponds at the 
theory of ^-bases for rational curves and our definition is equivalent to the usual 
definition as in [Chen, Wang(2003)[ Theorem 3, Condition 3]. 



However, the approach used here is actually inverse to the approach in the 
cited papers. In the latter the ruled surface is defined by a parametrization and 
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then the module of moving planes is studied, whereas here we fix a module that 
"looks like" such a moving plane module and then study the (generalized) ruled 
surface that corresponds to it. Note that by definition of the subspaces L(M,t) 
any element C of M can be considered a moving plane following Sm > m the 
sense that for all x € Sm there is a parameter t such that C(t) ■ x = 0. 

Note that A A B defines the so-called Plucker curve V in p<*(<*-i)/2-i ^ y 

p v : P 1 -~ ► pd(d-i)/2-i 

* ^ ([1,2] : [1,3] : ... : [d-l,d]) 

where [i,j] = AiBj — AjBi. We will denote k = deg ip-p the degree of the 
parametrization, which is the cardinality of the fiber of a generic point in the 
image of pp . Note that (p-p and k are the same for any choice of quasi-generators 
of M. 

Proposition 4. For any pair of quasi-generators A, B of M we have the degree 
formula 

k ■ deg S M = deg(A A B) — deg q A ,B, 
where qa,b — gcd(AA-B) and k — degpp. Moreover, we have degV — deg S 



M ■ 



Proof. The proposition and the proof are similar to Lemma 1 in Chen et al.(2001)] 
and to Theorem 5.3 in Pottmann et al.(1998)| . 



The implicit degree of the hypersurface Sa,b is the number of intersections 
between a generic line and the hypersurface. The generic line L(s) is defined by 
two points in the space L(s) = Hq + sHi, where Hi = (ha, ha, hid), i = 0, 1. 
The line L(s) intersects the hyperplane {A(i) ■ x} if and only if Hq ■ A(t) + sH\ ■ 
A(t) = 0. Since the line L(s) should intersect the hyperplane {B(t) ■ x} too, we 
see that the implicit degree is the number of intersections of two curves in the 
(t, s) plane: 

H ■ A(t) + sHx ■ A(t) = 0, 
H ■ B(t) + sHi ■ B(t) = 0. 
Eliminating s from the above equation we have 



H -A(t) H x -A{t) 
Ho -B{t) HfB(t) 



(Hq A H±) ■ (A(t) A B(t)) = 0, (4) 



where C ■ D means a standard scalar product of two vectors C,D S R d ( d_1 )/ 2 . 
The number of solutions of ([4]) is the number of intersection points of the Plucker 
curve with a generic hyperplane in p <J ( <i - 1 )/ 2 - 1 ; so degV = deg<S>M- 
Now it is known (see for example |Dohm(2006) Theorem 1] ) that 



k ■ deg V = deg(A A B) - deg qA,B 

and the proposition follows. □ 

We have yet to show the existence of the /i-basis. To this end, we propose 
an algorithm for its computation, the basic idea of which is to reduce qA,B — 
gcd(A A B) to a constant using the so-called Smith form of the 2 x d matrix 



W A . 



A x A 2 ■■■ A d 
Bi B 2 ■ ■ ■ Bd 
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and then render the leading vectors linearly independent by a simple degree 
reduction. The Smith form is a decomposition Wa,b — U -S-V , with unimodular 

U e R[i] 2x2 , V € R[t] dxd , and 

s =(l o ^ °^ntr d 

\0 q A ,B ••• OJ 

It always exists and can be computed efficiently by standard computer algebra 
systems. 

Algorithm (/i-basis) 

1. INPUT: Quasi-generators A = (Ai, A 2 , . . . , A d ),B = {B x , B 2 , B d ) E 
R[t] d of the module M 



2. Set 

W A . B 

3. Compute a Smith form 



A x A 2 ••• A d 
B\ B 2 ■ ■ ■ B d 



w AB = u-( 1 •'• °Vv 

A ' B \0 q A ,B •■■ 0^ 

with unimodular U £ R[t] 2x2 ,V 6 M[t] dxrf . 

4. Set W to be the 2 x d-submatrix consisting of the first two rows of V. 

5. If the vector of leading terms (with respect to the variable t) of the first row 
is h times the one of the second row, h £ M[t], set W :— i\ J 1 ] ■ W . 



v° 1 , 

6. If the vector of leading terms (with respect to the variable t) of the second 
row is h times the one of the first row, h G M[t], set W' := ( ^ ^ 

7. If the preceding two steps changed W go back to Step 5. 

8. Set A, B to be the rows of W. 

9. OUTPUT: A /i-basis A, B of the module M 

As we shall see in Section[5j the case we are interested in is the case d — 6, so we 
are dealing with very small matrices and the computations are extremely fast. 
Note that we actually only need the first two rows of V, so we could optimize 
the algorithm by modifying the Smith form algorithm used as not to compute 
the unnecessary entries of the matrices U and V. Generally, the number of 
elementary matrix operations in Step 5 and 6 is very low. In the worst case, it 
is bounded by the maximal degree of the entries of the matrix W in Step 4 of 
the algorithm, since each step reduces the maximal degree in one of the rows of 
W. 

Next, we will show that the output of the above algorithm is a /i-basis and 
that the resultant of a /i-basis A, B of the module M — (A, B) is an implicit 
equation of Sm ■ In Section [5J we will use these results for a special choice of A 
and B to compute the implicit equation of a canal surface. 
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Lemma 5. The output of the above algorithm is a fi-basis and we have k ■ 
dcgiSj\/ = deg M, where k = degip-p. 

Proof. Let A(t), B(t) be the output of the above algorithm. By construction it is 
clear that A(t), B(t) are quasi-generators of M and that q~A,B = gcd(A A B) = 1. 
Furthermore, we have deg(A A B) — deg(A) + deg(B), because the vectors of 
leading terms of A(t) and B(t) are linearly independent. So by Proposition [3] 
we deduce 

k ■ deg S M = deg(i A B) - deg(g) 
= deg{A A B) 
= deg(A) + deg(B) 

Moreover, by definition we have degM < deg(A A B) and if A, B are quasi- 
generators such that deg(A A B) is minimal, the degree formula gives deg(A A 
B) = dcg(yl AS) — deg(g) < deg M, which shows that k ■ deg Sm = deg M, and 
as a consequence that A(t), B(t) is indeed a /i-basis. □ 

Proposition 6. Let A(t),B(t) be a fi-basis of M and x = (x±, x%, . . . , Xd) T 
variables. Then 

Res t (A(t)-x,B(t)-x) = Fl M 
where F$ M is the implicit equation of the hypersurface Sm ■ 

Proof. First, we will show in the same way as in [Dohm(2006)| Theorem 9] that 
Res t (A(t)-x,B(t)-x) is geometrically irreducible, i.e. the power of an irreducible 
polynomial. As we shall see in Proposition the intersection of the hyperplancs 
{A(t) ■ x} and {B(t) ■ x} is of codimension 2 for any parameter t £ P 1 . So the 
incidence variety 

W = {(t,x) e P 1 x P^liO) • x = B(t) ■ x = 0} C P 1 x P^ 1 

is a vector bundle over P 1 and hence irreducible. So the projection on P d_1 
is irreducible as well and its equation, which is by definition the hypersurface 
defined by Rest(A(t) ■ x, B(t) ■ x), is a power of an irreducible polynomial. 

As we have remarked earlier, the resultant of two quasi-generators is always 
a multiple of the implicit equation of Sm, so Rest(A(i) • x, B{t) • x) is a power 

But using the degree property above we see 

deg(Res t (i(t) • x,B(t) ■ x)) = deg(i) + deg(B) = k ■ degS M 
which implies that Rest(>l(t) ■ x, B(t) ■ x) equals ^5 ■ □ 

Remark 7. It is known that the Pliicker curve V can be properly reparametrized, 
i.e. there exists a rational function h of degree k such that A A B = C o h, 
where C is a proper parametrization of V . It is tempting to use this proper 
reparametrization in order to represent the implicit equation F$ M of Sm di- 
rectly as a resultant as in the proof of |Dohm(2006)| Theorem 3]. However, h 
does not necessarily factorize A and B, i.e. it is not sure that there exist A' and 
B' with A = A' o h and B = B' o h, which would be needed to do this. 
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In the following we present some properties of /x-bases. Note that the prop- 
erties in Propositions [SJ H] are similar to Chen, Wang(2 003)| Theorems 1,3. 



However, we give different proofs by deducing them from the degree formula 
and Lemma [5] 

Proposition 8. Let M = (A, B) and let A,B be a fi-basis of the module M. 
Then the following properties hold: 

1. The vectors LV(A), LV(B) are linearly independent. 

2. A(to),B(to) are linearly independent over C for any parameter value to £ C. 

Proof. 1. If LV(A), LV(B) were linearly dependent, this would imply that 
k ■ degSM = deg(A AB) - deg(q^ g) < deg(A) + deg(B) = degM which is a 
contradiction to Lemma 

2. Suppose that A(to),B(to) are linearly dependent for some to £ C. This 
is equivalent to saying that the matrix ^ is not of full rank, which means 
that all 2-minors vanish. So to is a root of ^ and as above we deduce k • 
degiSjvf = deg(A A B) — deg(q^ g) < deg(A) + deg(B) — degM which is again 
a contradiction to Lemma [5l □ 

Proposition 9. Let M = {A, B) and assume that A, B satisfy conditions 1,2 
from Proposition Then any element D £ M has the following expression: 
D = hi A + h2B for some hx,fi2 G R[i], i-e. A, B are generators of the module 
M over the polynomial ring R[t]. Moreover, the pair A, B is a fi-basis of the 
module M . 

Proof. Let D 6 M , it can be expressed as 

with a, b, c, d € R[f] and co-prime numerators and denominators in the rational 
functions ? and %. Furthermore, we may assume that gcd(a,c) = 1, because 
if gcd ^ c ) is a linear combination of A, B, then so is D. Multiplying both 

sides of the above equation with bd we obtain bdD = adA + bcB or equivalently 
b(dD-cB) — adA and since b divides neither a nor A (if it divided A, for any root 
to of b and any constant a we would deduce the relation = aA(t ) + • B(t ), 
which contradicts property 2 in Proposition [5]) , one concludes that b divides d 
and by a symmetric argument that d divides b, so we may assume b = d. So we 
have 

bD = aA + cB 

and plugging a root to of b into the equation, we would obtain a non-trivial linear 
relation between A and B, again a contradiction to Proposition [5] This implies 
that b and d are constant, which shows that any D € M can be expressed as 
linear combination of A and B over M.[t}. In other words: A and B are not only 
quasi-generators of M, but actually generators in the usual sense, i.e. over R[t}. 

Suppose that dcgA < dcgB and let M = (Pi,P2). Then we proved that 
Pi = hiiA+h i2 B, i = 1, 2 for some polynomials hij £ R[t\. Since LV(A), LV{B) 
are linearly independent LV(hnA) and LV(hi2B),i = 1,2 do not cancel each 
other. Therefore, degPi > degB (if h u j^ 0) or degP 2 > degi? (if h 22 + 0). 
Also degPi > deg A and degP 2 > deg A So, we see that degPi + degP 2 > 
deg A + deg B, i.e. a pair A, B is a /x-basis of the module M. □ 
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3 Elements of Lie and Laguerre Sphere Geome- 
try 

Here we shortly recall the elements of Lie and Laguerre Sphere Geometry (cf. 
[Cecil(1992)[|Pottmann, Peternell(1998)||Krasauskas, Maurer(2000)| ). We start 
from the construction of Lie's geometry of oriented spheres and planes in R 3 . 
Let p 6 R 3 , r £ R. The oriented sphere S p<r in R 3 is the set 

S P , r = {vel 3 |(v-p)-(v-p)=r 2 }, 

where by v-w we denote the standard positive definite scalar product in K 3 . The 
orientation is determined by the sign of r: the normals are pointing outwards 
if r > 0. If r = then S P fi = {p} is a point. Let neR 3 with n • n = 1 and 
The oriented plane Pn : h in R 3 is the set 

Pn,h = {v £ M 3 |v- n = h}. 

The Lie scalar product with signature (4, 2) in R 6 is defined by the formula 

— X\Z2 — 

[X, Z\ = h + + X5Z5 - XqZq. 

for x = (xi, . . . , Xq) and z — (z\, ... , zq). In matrix notation we have 

[x,z]=xCz T : where xC = (—0:2/2, —x\/2, X3, X4, X5, —xq). (5) 

Denote y — (u : yo : y\ : yi : j/3 : t/4) £ P(R 6 ) = P 5 and define the quadric 

Q = {y £ P 5 I [y,y] = -uy Q + yf+ y\ + y\ - y\ = 0} (6) 

where [., .] is the obvious extension of the Lie scalar product to P 5 . Q is called 
Lie quadric. 

We represent an oriented sphere S Ptr (or an oriented plane P n ,h) as a point 
Lie(S Pi r) (resp. Lie(P nt h)) on the Lie quadric: 

Lie{S^ r ) = (2(p • p - r 2 ),2,2p, 2r) £ Q, peR 3 ,reM, 
Lie(P nJl ) = (2/i,0,n, 1) £ Q, n 6 M 3 , h £ R. 

It is easy to see that we have determined a bijective correspondence between 
the set of points on the Lie quadric Q and the set of all oriented spheres/planes 
in M 3 . Here we assume that a point q = (1 : : : : : 0) £ Q on the Lie 
quadric Q corresponds to an infinity, i.e. to a point in the compactification of 
R 3 . We say that q = (1 : : : : : 0) is the improper point on the Lie 
quadric. Notice that oriented planes in R 3 correspond to points Q f] T q , where 
T q = {V = ( u '■ Vo ■ 2/i : 2/2 : 2/3 : 2/4) G P 5 \y = 0} is a tangent hyperplane to the 
Lie quadric at the improper point q. 

Two oriented spheres S Pliri , S P2tT2 are in oriented contact if they are tangent 
and have the same orientation at the point of contact. The analytic condition 
for oriented contact is 

llPi -P2II = ki -r 2 |, 
where ||p 1 — p> 2 1| denotes the usual distance between two points in the Euclidean 
space R 3 . One can check directly that the analytical condition of oriented 
contact on the Lie quadric is equivalent to the equation 

[Lie(S Pl , ri ), Lie(S P2tr2 )} =0. 
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It is known that the Lie quadric contains projective lines but no linear sub- 
spaces of higher dimension (Chapter 1, Corollary 5.2 in Cecil(1992)| ). More- 



over, the line in P 5 determined by two points k\, k 2 of Q lies on Q if and only 
[fci, k 2 ] = 0, i.e. the corresponding spheres to k\, k 2 are in an oriented contact 
(Chapter 1, Theorem 1.5.4 in [Cecil(1992)] ). The points on a line on Q form so 



called parabolic pencil of spheres. All spheres which correspond to a line on Q 
are precisely the set of all spheres in an oriented contact. 

Remark 10. Here we use a slightly different coordinate system in Lie Geometry 
than in the book [Cecil(1992)| . The scalar product as in [Cecil(1992)| may be 
obtained applying the following transformation: 

x[ = (Xl + X 2 )/2,x' 2 = (x 2 - X 1 )/2,x' 3 = £3, £4 = ^4,^5 = X 5 ,x' 6 = X 6 . 

We show now that the set of points y in Q with ?/o ^ is naturally diffco- 
morphic to the affine space Mr. This diffcomorphism is defined by the map 

0: Q\T q -» M 4 , 

(u : y : Vl : y 2 : y 3 : Vi ) ~ (H^^), 

where T q = {y = (it : y : yi : y 2 : 2/3 : Vi) e P 5 | y = 0} as before, i.e. the 
tangent hyperplane to the Lie quadric Q at the improper point q = (1 : : : 
0:0:0). Let v — (vi, v 2 , v 3 , V4), w — (wi,w 2 ,W3, W4) £ P 4 and denote by 

(v, w) = V1W1 + V 2 W 2 + V 3 W 3 - V4W4 

the Lorentz scalar product on P 4 , which can be seen as the restriction of the Lie 
scalar product [., .] to R 4 . The affine space R 4 with the Lorentz scalar product 
is called the Lorentz space and denoted by Rf . 

Let y = (2/1,2/2,2/3,2/4) £ R ■ One can check that inverse map of 4> is given 
by the formula: 

r 1 (2/) = ((2/,2/),i,2/)eQ\r, 

Notice, that 4>(Lie(S Ptr )) = (p,r), i.e. the sphere S Pyr 6 R 3 corresponds to a 
point (p,r) G R 4 . The map <j) can be extended to a linear projection $ from 
Q \ {q} to P 4 denned as 

$: Q\{q} -» P 4 

(m : 2/o : 2/1 : 2/2 : 2/3 : 2/4) (2/0 : 2/1 = 2/2 : 2/3 : 2/4) 

The points of QC\T q can be represented as Lie(P n j l ) = (2h, 0, n, 1) and 
these points correspond to planes in R 3 . Note that 

$(Lie(P n , h )) = (0, n, 1) G n = {2/0 = 0, y\ + 2/ 2 2 + 2/ 3 2 - vl = 0} 

are infinite points to the natural extension of R 4 to P 4 which correspond to a 
pencil of parallel planes in M 3 . The quadric £1 is called absolute quadric. The 
preimage of the map $ has the following form 

®~ 1 (V) = ((2/, y) ■ 2/0 : voyi ■ 2/02/2 : 2/02/3 : 2/02/4) (7) 

where y = (y : y x : y 2 : 2/3 : 2/4) G P 4 and y = (2/1,2/2,2/3, 2/4) as before. 
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A direct computation shows that for v,w G R 4 

-2[cj)- 1 (v),(j)- 1 {w)] = (v-w,v-w) (8) 

The formula shows that two oriented spheres defined by v, w (i.e. spheres 
S/ Vl ,v2,va),v4 an d S(tm,i>i2,wa),w4.) are m oriented contact if and only if (v — w, v — 

w) = 6. 

Let us define two maps: an embedding id : R 3 — > R 4 ,id(p) = (p, d),d E R 
and a projection it : R 4 — > R 3 ,7r(p, r) = p, where r e R. We will treat points 
io(R 3 ) as spheres with zero radius and identify them with R 3 . All interrelations 
between the spaces introduced above can be described in the following diagram 

Q\T q C Q\{q} C P 5 

R 3 % R 4 c P 4 (9) 

II I* 

R 3 = R 3 

Definition 11. For an oriented surface (curve or point) M C R 3 define an 
isotropic hypersurface Q(A4) C P 4 as the union of all points in R 4 which cor- 
respond to oriented tangent spheres of M. Let Qd(M) — Q{M) (~l {2/4 = dyo} 
be a variety which corresponds to tangent spheres with radius d of M.. The 
set Env ( i(A / () = ^{Qd{M.)\^) C R 3 are centers of spheres with radius d tan- 
gent to M. The set Envd(.M) is called d-envelope of the variety M. Since 
0(M) = \J d Qd{M) we can treat the isotropic hypersurface Q(M) as the union 
of all d-envelope to the variety M.. 

If y, a £ R , ao, yo G K, we define a function 

g((a : a), (y : y)) = (ay - a y, ay - a y) = 2/o a o ( — -—,— - — (lO) 

\ao yo «o yo/ 

= Va(a,a) -2a Q y (a,y) +al(y,y). 

Let (j/o : y) be such that g((ao ■ a), (yo : y)) = 0. By the formula ([5]) we see that 

spheres SV^ 03 £3 \ ^4 an d SV3/1 ^ jj\ M are in oriented contact. Therefore, in 

\ "0 ' "0 ' "0 / ' "0 v yo ' vo ' wo ) ' va 

the same manner as previously, we define the isotropic hypersurface G((a,o : a)) 

as follows 

g((a : a)) = {(y ,y) G P 4 | g((a : a), (y : j/)) = 0} C P 4 . 

In fact, Q((oq : a)) is a quadratic cone with a singular point at a vertex (ao : 
a) G P 4 and may be viewed as the set of all spheres which touches the fixed 

sphere S / ai a 2 a 3 \ ai ■ After the restriction to the linear subspace 2/4 = dj/o this 

V a o ' a o ' a o J ' "-a 

hypersurface consists of all spheres with radius d which are in oriented contact 

with the sphere S/ ai a2 a 3 \ a 4 which we denote as Gd((&o ■ a)) = ^((o-o : 

V a o ' a o ' a o J ' a o 

a)) n {2/4 = dyo}- We notice that Od{{ao : a))\ Vo= i is defined by the equation 
(ai - a yi) 2 + (a 2 - a y 2 ) 2 + (a 3 - a y 3 ) 2 = (a 4 - a d) 2 , i.e. 



TT(Gd((a a : a))\ m 4) = g/ 01 „ 2 „ 3 n „ 4 _ ao „ = Env_ d I g/ ai a2 , 3 n « 4 ) and 

V a o' a o' a o/' "0 \ \°ll'"0'«0/' o 



a a a / 



</d((aa : a))| R 4 = i d I S 



11 



Therefore, in this case, the isotropic hypersurface Q((ao ■ a)) may be treated 
as a union all envelopes to the sphere S / ^ ^ ^4 • In the next section we 

generalize the definition of the isotropic hypersurface Q(M) for a curve M in 
R 4 (or P 4 ). 

All lines in Rf with directional vectors v can be classified into three types 
depending on the sign of (v,v): (+)-lines, (O)-lines (also called isotropic lines), 
and (— )-lines. 

4 The isotropic hypersurface and d-envelopes 

In this section, we will see that the definition of the canal surface is not obvious 
and and we will introduce some geometrical object related to it. A canal surface 
is given by a so-called spine curve £, which is the closed image (with respect to 
the Zariski topology) of a rational map 

R — > R 4 

1 ,_> ( e i(*) £l(t) es(t) e 4 (t) \ 
1 ^ Veo(t)' e (t)> e (t)> e (t) J 

with polynomials eo(i), . . . , e^t) £ R[t] such that n — maxi = o,..4{deg(ei(i))}. 
For abbreviation, we usually skip the variable t in the notations. The spine curve 
describes a family of spheres (5/ ei (t) e 2 w sw\ g^m | i G R} whose centers are 

given by the first three coordinates ( ei i t l , 62 ¥1 , ) and whose radii are given 

° J \eo(t) ' e (t) ' e (t) y to 

by the last coordinate f^jj- Intuitively, the canal surface is the envelope of this 
family of spheres, but there are some subtleties to consider before we can make 
a precise definition. 

We can also consider the spine curve as a projective curve £ given as the 
closed image of a parametrization 

P 1 —4 P 4 

t ^ (e (t) : ei(t) : e 2 (t) : e 3 (t) : e 4 (*)) 

with the non-restrictive condition gcd(eo, ■ ■ ■ ,e 4 ) = 1, which means that there 
are no base-points (i.e. parameters for which the map is not well-defined). 

Note that in this case the polynomials e, are actually to be considered as 
homogenized to the same degree n with respect to a new variable s. As there 
is a one-to-one correspondence between the univariate polynomials of a certain 
degree and their homogeneous counterparts, we will keep the notation from 
above and distinguish between the affine and projective case only where it is 
necessary to avoid confusion. 

In the following we use the notations 

e = (ei,e 2 ,e 3 ,e 4 ), y = (2/1,3/2,2/3,2/4), 

e = (e : ei : e 2 : e 3 : e 4 ), y = (y : 2/1 : 2/2 : 2/3 : 2/4 )• 

We first proceed to define a hypersurface in P 4 which is closely related to the 
canal surface. 
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Definition 12. The isotropic hypersurface Q{£) = {y | G(y) = 0} C P 4 asso- 
ciated with the (projective) spine curve £ is the variety in P 4 defined by the 
polynomial G(y) = Res t(.gi,92) where 

9i(y,t) = g{e,y) = {e yi - eij/ ) 2 + {eoV2 - e 2 yo) 2 + 
(e a y 3 - e 3 y ) 2 - (e y 4 - e 4 y Q ) 2 
= (eoV ~ Voe, e y - y Q e) = el(y,y) - 2(e e, y y) + yl(e, e), 

92{V,t) = dgi ^^ = 2 ( e oe' (y,y) - ((e e)', y y) + y$(e', e)). 

So, we define Q{£) as the envelope of the family of isotropic hypersurfaces Qie) = 
Q((eo(t) : e(t)). 



In the previous section we showed that Sti(e) |r4 = Env-^ ( S 1 
This interpretation leads to the following definition. 



Definition 13. The d-envelope associated with the (projective) spine curve £ 
is defined as the hypersurface Envd(£ ) C P 3 given by the implicit equation 

Gd{yo,yi,y2,y3) = Res t (gi\ yi =-d y „, g2\ Vi =-d Vo ) = ^s t (g 1 ,g 2 )\y i =-dy a , 

i.e. the equation obtained by replacing y 4 in G(y) by — dy , where d 6 K. 

The affine envelope Envd(£) at distance d is the restriction of Envd(£) to the 
affine space R 3 , defined by the equation Gd\ Vo =i = Res t (gi, g2)\ y4 =-dy ,yo=i> 1 - e - 
by setting y = 1. 

So G(£) contains all offsets associated with the spine curve £. Indeed, the 
surface 

Env d (£) = g(£) n {y 4 = -dy } 

is a hyperplane section of G(£), which can be interpreted as a parametrization 
of all offsets (with respect to the parameter y 4 ). 

The special case d = is particularly important. For the real part of Envo(£ ) 
to be non-empty, one has to suppose that £ has tangent (+)-lines almost every- 
where, or equivalcntly that (e, e) > almost everywhere. Env (£) is the en- 
velope of the family of spheres in R 3 given by the spine curve £ and Envd(£ ) 
is the envelope of the same family of spheres with radii augmented by d. For 
instance, circular cylinders or circular cones (call them just cones) are envelopes 
Env (£) of (+)-lines C and vice versa. In the literature, the canal surface C is 
usually defined as this envelope Env (£). However, we will show in an exam- 
ple that these envelopes can contain "unwanted" extraneous factors, which are 
geometrically counterintuitive. 

Example 14. Consider the spine curve £ given by 

' ei(t) e 2 (t) e 3 {t) e 4 {t) \ f l-t 2 2t 1 
e (t) ' eo(t) ' eo(t) ' e (t) ) \ 1 + t 2 ' 1 + t 2 ' ' 2 

The first three coordinates describe a circle in the plane and moving spheres of 
constant radius along this curve, so intuitively the envelope should be a torus 
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T. But it turns out that the implicit equation of Envo(£) is up to a constant 
computed as 

G = Rfis t (ffi,52)| w=0lWJ= i = (y? + 2/ 2 2 ) 2 (4y 2 + Ay 2 2 + Ayj + 8y± + 3)F r 

where Fq- is indeed the equation of the torus. To understand where the other 
factors come from, consider the following: For a given parameter t, the equations 
gi and gi define spheres S\ (t) and Si (t) in R 3 and 

Env (£) = U^i(t)n5i(t) 

t 

of the intersections of these spheres (actually this is nothing else than the geo- 
metric definition of the resultant). Now, while for almost all t this intersection 
is a transversal circle on the torus (often called characteristic circle in the lit- 
erature), it can happen that the spheres degenerate either to planes or to the 
whole space. In our example, for the parameters t = i and t = — i we have 
9i(i) = 9i(—i) = 0; 92(1) = —iyi + D2 and g2(—i) = iyi + 2/2, so the intersection 
in those parameters actually degenerates to (complex) planes which correspond 
to the factor (—iyi + 2/2) (%i + 2/2) = 9% + 2/f- In the parameter t = 00, both 
<7i and gi define the same sphere whose equation 4yf + 4y| + 4j/| + 8yi + 3 is 
the other extraneous factor. This kind of phenomenon can also happen for real 
parameter values, but it is interesting to remark that even though we consider 
a real parametrization, non-real parameters can interfere with the envelope, 
because the resultant "knows" about them. 

This example shows that Envo(£ ) is not a suitable definition for the canal 
surface C and we will later develop one that avoids the kind of extraneous 
components we have observed. 

Remark 15. Sometimes, in the literature, Envo(£) is defined in affine space as 
the resultant 

60(2/1,2/2,2/3) = Res * (9i (2/1 , 2/2, 2/3,*), h (2/1, 2/2, 2/3, *)), where 
9i = eg/i, 52 = eo/ 2 , 

\2 / x 2 / x 2 / \ 2 

ei \ , / e 2 \ , / e 3 \ / e 4 



h = 2/i + 2/2 +2/3 , 



h 



dh 
dt 

or in other words by deriving the affine equation of the sphere after the substi- 
tutions 2/4 = 0, 2/0 = 1 and homogenizing afterwards. Note that in this case fi 
and gi are linear in J/1,2/2,2/3- Let /1 = ffi| B4 =o, OT =i and /s = 52| w =o,y =i- An 
easy computation shows that we have the following equalities 

9i = hy h = eo/2 - 2e /i. 
Therefore, by properties of the resultant (|29l) . (|30|) . we have 

R.es t (<h,(/2) = Res t (/i,e / 2 - 2e /i) 
= Res eo/2) 
= Rcs t (/i,e ) • Res t (/i,/ 2 ). 

Hence, we have Go = Res t(fi, eo) -Go, so there are even more extraneous factors 
than before due to the roots of 6q. 
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Linearizing the problem 

The main idea to understand and eliminate the extraneous components that 
appeared in the example is to linearize the equations g± and g 2 by replacing 
the quadratic term (y, y) by a new variable u (or more precisely uyo to keep 
the equations homogeneous). This will make the results developed in Section 2 
applicable. Geometrically, this means that we will pull back the spine curve to 
Q via the correspondence 4>. 

For a spine curve £ E M.f we define a proper pre-image £ in the Lie quadric 
Q as the closure of the set £ = <J> _1 (£) in Q. It is immediate by ([7|) that the 
parametrization of £ is 

pi —4 Q c P 5 

t h-> ((e(t),e(t)) : eg(t) : e (t) ei {t) : e (t)e a (t) : e (t)e 3 (t) : e (i)e 4 (t)) 

(11) 

We can now define the envelopes associated with this new spine curve as follows. 

Definition 16. The variety TL(£) C P 5 associated with £ is the hypersurface 
in P 5 defined by the implicit equation H(y) = Rest(/ii, h 2 ) where y = (u : yo ■ 
Vi ■ 2/2 ■ 2/3 : 2/4) and 

hi(y,i) = -2[y,£(t)]=uel+y (e,e) -2{e e,y), 

h 2 (y,t) = dhl ^ t] = -2[y,£'(t)} = 2(ue e' +y (e',e) - ((e e)',y)). 

Similarly, the variety TLd{£) C P 4 is defined by the implicit equation 

Hd(u,y Q ,yi,y 2 ,y3) = Res t (/ii| y4= _ dW) , h 2 \ yi =-dy ) = Res f (hi, h 2 )\ yi =-dy , 

i.e. the equation obtained by replacing j/4 in H(y) by —dyo, where d £ M. 

Of course this is nothing else than substituting (y, y) in g\ and g 2 by uyo and 
dividing by y , so = h l ((y,y),y^,y y),i = 1,2, i.e. & = hi o $ _:L ,i = 1,2. 
Now as an immediate corollary we obtain 

Proposition 17. FF^/i the notations as above we have 

G(y) = H((y,y),yly y), i.e. G = H o (12) 

and 

Gd{yo,yi,y2,ys) = H d {yl + yl + yl~ d 2 yQ,yo,y y 1 ,yoy 2 ,yoy3)- (13) 

To sum up, we have defined two hypersurfaces as resultants of two quadratic 
forms: Env^) C P 3 , which are the offsets to the spine curve £ , and Q{£) C P 4 , 
which can be interpreted as a parametrization of those offsets. As seen in an 
example, these definitions can lead to additional components which are against 
the geometric intuition, so it is desirable to give another definition which avoids 
those extra factors. To this end, we have linearized the problem by replacing the 
quadratic polynomials g\ and g 2 by linear forms h\ and h 2 by substituting the 
quadratic term by a new variable and have seen how to reverse this substitution. 
Geometrically, this means that we replace the hypersurfaces Envd(£) and G{£) 
by hypersurfaces 7id(£) and 7i{£ ) in one dimension higher. 

This has the advantage that we can now apply the technique of /i-bases 
developed earlier to understand and eliminate the extraneous factors of TL d {£) 
and H(£) and then come back to P 3 (resp. P 4 ) with the substitution formulae 
of Proposition [T71 
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5 The dual variety, offsets, and the canal sur- 
face. 



In this section, we will finally be able to define the canal surface C (and more 
general offsets to it) and the so-called dual variety r(£), which can be seen as 
a parametrization of the offsets to C. 

Up to the constant —2 the system hi = hi = is equal to 

'a A it) 

where the matrix C is defined by the formula ([5|). 

We can interpret the variety Ti.(£) defined by (114|) as a dual variety to the 
curve £ with respect to the Lie quadric Q, i.e. the dual variety to the curve 
£(t)C. Indeed, this dual variety consists of the hyperplanes which touch the 
curve £{t)C. The first equation in (|14[) means that the hyperplane contains 
the point £(t)C, the second equation means that the hyperplane contains the 
tangent vector £'(t)C to the curve £{t)C. 

In order to simplify notation we denote 



£{t)Cy T = 0, 
£'(t)Cy T = 0, 



(14) 



E = £{t)C=\ — -, , eoei,eoe 2 ,eoe3,-eoe 4 J (15) 

E' = £'(t)C=(—eoe' ,—{e',e),eQe 1 + e e' 1 ,e' e2 + e e' 2 , 

e'o e 3 + e e 3 , -ege 4 - eoe^) (16) 

and we have that Tt(£) — Se,e' by ©. As we have seen in Section 2, this surface 
contains extraneous factors which correspond to the roots of the 2-minors of the 
matrix We,e' , but which can be eliminated by replacing E, E' by a /z-basis of 
the module (E, E'). It is thus natural to make the following definition. 

Definition 18. We define the dual variety V{£ ) C P 5 to the curve £ as the 
hypersurface 

V{£) = S {E , EI) (17) 
where (E, E') is the module quasi-generated by E and E' . 

By the results of Section [H it is immediate that V(£) C Ti{£) does not 
contain the components of 7i(£) caused by parameters t where WE,E'(t) is not 
of full rank or equivalently, where the intersection of the hyperplanes defined by 
hi and h 2 is of codimension 1, i.e. the hyperplanes coincide. So we can deduce 

Proposition 19. Let Ei,E 2 be a fi-basis of the module quasi- generated by E 
and E' and let k be the degree of the parametrization E A E' as in Section [Jl 
Then 

k ■ deg V{£) = deg E x + deg E 2 = deg(E A E') — deg q E ,E> , 
where qE,E' = gcd(E A E') and 



Res t (Ei-y T ,E 2 -y T )=F* {£): 



where i 7 ", 



v(£) * s ^ e implicit equation ofV(£). 



16 



Proof. It follows directly from Propositions 0] and [5] □ 

Of course, the same considerations can be applied to the hypersurfaces TLd{£) 
and we make the analogous definitions. Substituting 2/4 = — dyo in hi and h 2 
corresponds to replacing E and E' by two linear forms 

( e 2 (e, e) \ 
D = ( —|, '- \- deoe 4 ,eoe 1: e e2,e e 3 J (18) 

£>' = (-e e , -(e', e) + d(e e 4 - e e' 4 ) : e' ei + e ei, e' e 2 + e e' 2 ,e' e 3 + e e' 3 ) 

with D,D' G R 5 . Now 7id(£) — Sd.d 1 and one makes an analogous definition: 

Definition 20. We define the hypersurface Vd{£ ) as 

V d {£) = S {D , D , } c P 4 (19) 

where {D,D'} is the module quasi-generated by D and D 1 . 

In this case also, Vd(£) C TLd{£) does not contain extraneous factors due to 
the parameters t where the rank of Wo.D'it) drops. At this point, it should be 
remarked that while we clearly always have 

V d (£) C V(£) n {y 4 = -dy } 

this inclusion is not necessarily an equality (note that we had Env^f) = Q(£) n 
{i/4 = — dyo} for the corresponding varieties). Analogously to Proposition [T9l 
the following holds. 

Proposition 21. Let D\,D2 be a [i-basis of the module quasi- generated by D 
and D' and let k be the degree of the parametrization DAD' as in Section [Jl 
Then 

deg Vd{£) = deg D x + deg D 2 = deg(L> A D') - deg q D>D , , 
where qD,D' — gcd(D A D') and 

Res t (L>i • (u,y ,y 1 ,y2,y3) T ,D 2 ■ (u,y ,y 1 ,y2,y 3 ) T ) = F v d (£) 
where F Vd ^ is the implicit equation ofVd(£)- 

Proof. It follows directly from Propositions 0] and El □ 

Finally, we can use the correspondance of Proposition [T71 to define the canal 
surface. 

Definition 22. The T-hypersurface is defined as 

r(£) - $(v(£) n Q), 

and the offset OSd(£ ) at distance d to the canal surface C is 

OK d (£) = <i> (V d (£)nQ d ), 

where Q d = {(u : y : y x : y 2 : y 3 ) £ P 4 - uy a + y\+yl + yl~ d 2 yl = 0} and 
^(w, yo: 2/1, 2/2, 2/3) = (2/0,2/1,2/2,2/3)- The canal surface itself is the special case 
d = or in other words C — Offo(£). 
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Note that the extraneous factors of H d (£) and Ti{£) are in one-to-one cor- 
respondence with the extraneous factors of the corresponding hypersurfaces 
Env f j(£) and Q{£) since they are caused by parameter values where the in- 
tersection of h\ and /12 (resp. g\ and (72) is of codimension one. So T(£) and 
Cd{£) contain no such factors. 

In this section and the previous one, many different geometric objects have 
been defined. We illustrate in the following diagram how they are related in 
order to make the situation clearer. 



p4 p5 p5 p4 

u u u u 

g(£) *±- h{£)(iq 2 v(£)ng -±> r(£) 

U U U U (20) 

Env d (£) H d (£)f)Q d D V d {£) n Q d ^ OS d (£) 

n n n n 

p3 p4 p4 p3 



Note that the hypersurfaces in the third row are included in the correspond- 
ing hypersurfaces in the second row. The first column is the naive definition 
of the objects to be studied: F,nv d (£) is more or less a <i-offset to the canal 
offsets and Q{£) a hypersurface in one dimension higher containing all those 
offsets. However, they contain extraneous factors. So by passing to the second 
column, we linearize the hypersurfaces (i.e. we express them as resultants of 
linear forms) and can apply /it-bases to eliminate the extraneous factor, which 
gives the third column and finally go back down in dimension (by intersecting 
with Q and applying $ to obtain the objects we are interested in: the offsets 
Off d (£) (in particular the canal surface C = Offo(£)) and the T-hypersurface. 

5.1 The implicit equation. 

We can now describe how to compute powers of the implicit equations of the dual 
varieties V{£) and V d (£), the hypersurface T(£) and the offsets surface C d {£)- 
We should remark that these powers (which are the degrees of the parametriza- 
tions of the corresponding Pliicker curves) are in a way inherent to the geometry 
of the problem, as we shall illustrate in Example [271 They can be interpreted 
as the number of times the surface is traced by the spine curve. Note also that 
this not necessarily due to the non-properness of the spine curve: Even for a 
proper spine curve it can happen that the canal surface (or its offsets) is multi- 
ply traced, as in Example [771 

Algorithm (implicit equations) 

1. INPUT: A rational vector e(t) G R(t) 4 as in formula ©. 

2. Define E,E' G R[t} 6 as in formula (I5J) and D,D' G U[t} 5 as in formula 

3. Compute a /i-basis E\,E2 of the module (E,E') and a /i-basis Di,i?2 of 
the module (D,D') using the algorithm in Section [H 
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4. SetF y( g ) = Res t {E v y T ,E 2 -y T ) and F Vd ^ = Res t (D v (u, y Q , yi, y 2 , ys) T , 
D 2 ■ {u,y ,y 1 ,y 2 ,y 3 ) T ) = 0. 

5. Let F r(I) (y , 2/i, 2/2,2/3, y 4 ) = yoFy^^yl+yl+yj-ylj/yo, yo, 2/1,2/2, 2/3, 2/4), 
where k is a minimal integer such that -Fjvf) ^ s a polynomial. Similarly, 
set 

F C d (£) (2/0,2/1, 2/2, 2/3) = yo F v d (£)((yi +V2 +vl - d 2 y^)/y ,y ,y 1 ,y 2 ,y 3 ). 

6. OUTPUT: F v ,gy F v ,gy F r(s)' an< ^ F c d (£)^ wmcn are powers of the im- 
plicit equation of the varieties V(£), Vd(£), r(£) and C d {£) 

Note that the affine parts of these equations can be obtained by replacing 2/0 = 1 
before the resultant computation. 

5.2 The parametrization of the dual variety. 

We can describe the parametrization of V{£ ). The hyperplane defined by the 
equation 

det(y, € {t)C, £'{t)C, ai ,a 2 , a 3 ) = A lU + A 2 y Q + A 3Vl + ... + A 6 y 4 = (21) 

is tangent to the curve £(t)C , e R 6 ,i = 1,2,3 are three arbitrary points. 
By the definition a point on the dual variety V(£) is (Ai, ...,Aq), Define D = 
(£(t)C,£'(t)C, a%, 0,2, a 3 ) to be the 5 x 6 matrix with five rows £(t)C,£' (t)C, a\,a 2 , 
03. And let Di, i — 1, 6 be 5 x 5 matrices obtained from D by removing the 
i-th column. Then using the Laplacian expansion by minors for the first row of 
the determinant (|21[) we obtain the parametrization of V(£) as follows: 

c(D) = {det D l7 - det D 2 , det D 3 , - det D 4 , det D 5 , - det D 6 )/m C V{£), (22) 

where m — gcd(£*i, Dq). Here t, a\, a 2 , a 3 are arbitrary parameters. 

6 The implicit degree of the hypersurface 

The aim of this section is to get some formula for the implicit degree of the 
hypersurface T{£ ) in terms of the rational spine curve £ = {e(t) € P 4 }. Notice 
that the implicit degree of the canal surface C is less or equal than degT(£) 
because we always have the inclusion 

Off d (£) c T(£) n {2/4 = -dyo}, i.e. degOff d (£) < degT(£). 

So this formula gives upper bound for the degree of the canal surface. In the 
case of a polynomial spine curve the upper bound was obtained in the paper 
Xu et al. (2006)| . Note that for the computation of the implicit degree we do 
not need the implicit equation of the hypersurface. We believe that this formula 
is useful for higher degree spine curves because the computation of the implicit 
equation may be very difficult in practice. 

Let us remind that the pre-image $ -1 (r(£)) C Q is defined by the inter- 
section of two varieties V(£) fl Q. Let denote by G(u, yo, 2/1,2/2,2/3, 2/4) = the 
equation of V(£). The Lie quadric has the equation uyo = (y,y) (recall that 
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(y,y) — y\ +y\ + y| — y|)- By the definition (f2"2"|) the equation of is obtained 
after the elimination of the variable u from the equations of Q and V(£), i.e. 

r (£) : |^(2/o, 2/1, 2/2, 2/3,2/4) = VqG f^^> 2/0, 2/1, 2/2, 2/3, 2/4^ = oj , (23) 

where fc is a minimal integer such that the left side of the equation (|23l) is 
polynomial. We introduce the following weighted degree 

^(«' Cl 2/o fc2 2/^ 3 2/2 fe4 2/3 5 2/4 6 ) = 2k 1 + k 3 + k 4 + k 5 + k 6 , (24) 
d w (G(u, yo,2/i, 2/2, 2/3, 2/4)) = max{d u ,(m i )} 

where G = ^ c * m j * s a linear combination of the monomials rrij. Using this 
notation we have degT(£) = d w (G). 

Let us assume that the curve £ = {e(t, s) £ I 4 (t, s) S P 1 } has a homoge- 
neous parametrization. We introduce the following notations 

w = (wi,W2,w 3 ,W4,), Wj = e'jeo - eje' , j = 1,2,3,4, (25) 

7 = max{deg(w i ,t)}, (26) 

j 

where e\ means derivative with respect to t. We will say that the curve £ £ P 4 
is of general type if : 

gcd(u>i, u> 2 , w 3 , W4) = gcd(e ,eo) = gcd(e , (e,e)) = 1, deg£ = deg(e ,t) and 
the parametrization degree of the Pliicker curve ip-p : t — » E A E' is one, (27) 

(i.e. degcp-p = 1), where e = (ex, e2, e3, e 4 ) and i? as in formula (fT5"|) . 
The first equation in the system IT4^) has the following form: 

hi = £Cf = Ey T = -e>/2 - (e, e)y /2 + e (e, y). (28) 

The polynomial /ii is linear in the variables u, yo, 2/1, 2/2, 2/3, 2/4 an d has degree 2n 
in the variable t. The elimination of the variable t from the system h\ = h'i = 
is the reducible polynomial Rest(/ii, /i'jJ = Hi...HkG. By definition one of those 
factors is the equation of the dual variety V{£) — {G}. 

Proposition 23. If £ is a curve of general type then deg V{£ ) = An — 2, where 
n = deg£. Moreover, we have G ■ LC{h\) — Rest(/ii, h'i), where LC{h\) is 
the leading coefficient of the polynomial h\ with respect to the variable t and 
{y e P 5 I G{y) = 0} = V{£). 

Proof. For the curve of general type by the Proposition [19] we have 
degV(£) = degE A E' — deg(gcd(i? A E')). We can compute components of 
the Pliicker vector E A E' = ([1,2] : [1,3] : ... : [5,6]) G P 14 . For example, 
[1,2+.?] = eluij/2, j = 1,2, 3, 4 and [2,3] = (-(e, e){e ex)' + e ei(e, e)')/2. By 
assumption (127p we see 

gcd([l,3], [1,4], [1,5], [1,6]) = e 2 andgcd(e , [2,3]) = 1. Therefore gcd{EhE') = 
1. So we have deg V(£ ) = degi? A E' — An - 2. Notice that deg Rest {hi, h' x ) = 
An— I. Therefore we see deg G = deg Res* (/ii, h'-^/LCQii) = An - 2, i.e. G = 
is the implicit equation of V{£). □ 
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Thereinafter, we will show that for the curve of general type, we have 
d w (G) = 6n — 4, where n — deg£ . For this we consider another resultant 
Rest(/ii, /12) and show that d w (Rjes t (hi, /12)) = d w (G). We define h% in the fol- 
lowing way. Let g — e\ then we have the following equality h[g — h\g' = eohi, 
where 

hi = (2(e,e)eo-(e,e)'eo)yo/2+e (e'eo-eeo,y} = (2(e, e)e -(e, e)'e )2/o/2+e (w, y). 
In other words, /12 is the numerator of a rational function ( — 



We often use the following properties of the resultant 



Res(/i/ 2 ,/i) = Res(/i,/i)Res(/ 2 ,/i) (see [Cox et al.(1998)] , p. 73), (29) 



Rea(f,h) = cgr Res(/, r) (see [Cox et al.(1998)] , p. 70), (30) 



where h = qf + r, degr < deg ft- = m,f = a^x 1 + a\x 1 + ... + a;. 
We need an explicit formula for factors of the resultant 

LC^Restihuh^g - h ig ') =Res t (h u h[g), (31) 

where LC(hi) is a leading coefficient with respect to variable t of the polynomial 
hi. Indeed, since deg(/ii,i) = deg(g,t) then deg[h'^g — hig', t) + 1 = deg(/i' 1 g,t). 
Thus we obtain the formula in (|3~lj) from the property ([31?)) . The left side of the 
formula (|3"Tj) is equal to LC{h\) Res t (fti, eo/12) = LC{h\) Res(/ii, eo) Res t (/ii, /i 2 ). 
Otherwise, the right side of this formula is equal to Rest(/ii, /i^) Rest(/ii, eo) 2 = 
G- J LC(fti)Res t (fti,e ) 2 . 

Also, from the property J2U]) follows that Kest(h\, eo) = y^ 6 ^"'^. Therefore 
we have and d w (Rest(hi, /12)) = d w (G). In the lemma below, we prove that 
d w (Rest (fti, /i2)) = 6n - 4. 

We summarize our computations in the following 

Theorem 24. The degree of the hypersurface T(£) with the spine curve £ which 
satisfies the assumption \27\ ) is equal to 6n — 4 7 where n — deg£. 

Lemma 25. With the notation as above we suppose that the conditions are 
satisfied. Then we have the equality d w (Res t (hi, hi)) — 6n—4, where n — deg£. 



Proof. The weighted degree d w (G) may be viewed as a degree of a variety 
$({G} n Q), where $ : P 5 \ q — > P 4 is a linear projection from the improper 
point q — (1,0, ...,0) on the Lie quadric (see explicit formula ©). The degree 
of the variety ^({G} fl Q) can be computed constructively by counting points of 
intersection with a general line L £ P 4 . The pre-image of the line G := <f> _1 (iv) 
is a conic on the Lie quadric Q which passes through the improper point q. 
Hence the degree of $({G} n Q) is 2 deg G - i(q, {G} n G), where i(q, {G} n G) 
is the multiplicity of the intersection {G} fl G at the point q. 

We need a parametric representation of the general conic q G G C Q. As- 
sume that the conic G is in the parameterized plane P : q + k\u + fc2i>, where 
fci, fc 2 £ K 6 . The plane P intersects a singular cone {(a;, x) = x\+X2+x\—x\} on 
two lines. We choose two vectors k\, ki in these lines so that the first coordinate 
is zero, i.e. k% = (0, l,a),fc 2 = (0,1, b), a = (ax, a 2 , a 3 , a 4 ), b = (61, &2, ^3, 64) 
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such that (a, a) = (b, b) = 0. With the notations as above the general conic 
C := P fl Q has the following parametrization: 

C(u) := (ku — 1, ku 2 , ku 2 a + u(b — a)), where k = 2{a, b). (32) 

Since C(0) = q we can compute the multiplicity m = i(q, {Res(fti, ft,2, t)} PI C) 
as follows. Lets denote by cfti = h\\c, c/12 — h.2\c the restriction of polynomials 
hi, ft 2 to the conic C: 

c/ii = -el{ku-l)/2-kfu 2 /2 + e Q {e 1 ku 2 a + u(b-a)}, (33) 
cft 2 = u((2fe' — f'eo)ku/2 + eo(w, kua + b — a)), where /= (e, e). (34) 

The computation of the resultant gives Res t(chi, 0/12) = w" l (A + _Bu + ...), here 
m = {Res t (/ii, /12)} fl C). On the other side we can compute the number of 
common points (0, to) on curves {chi(u, t)} and {ch,2(u, t)} counted with multi- 
plicities. This number coincides with m (see |Buse et al.(2005)] , Proposition 5). 
The second curve cfi2(u,t) = u ■ ch^(u,t) is reducible. Therefore, the resultant 
with respect to t is 

Res(c/ii, ch 2 ) = Re$(chi,u)Res(chi,ch3) — u dcg( - hl ' t ' ) Res(chi,ch 3 ) 
= u 2n Res(chi,ch 3 ) 

The second factor has a representation ch 3 = uD(t) — N(t), where D(t) 
'If'] — f'eo)k/2 + eo(w, ka) and N(t) — eo(e'eo — ee' , a — b). It is easy to sec 
that gcd(N(t),el) = e . If cfti(0, t ) = ch 3 (0, t ) = then e (t ) = 0. The first 
curve {chi(u,t)} is hyper-elliptic, i.e. the projection to t axes pr : {chi} — > t 
is a map two-to-one. The curve {chi(u, t)} has the following discriminant with 
respect to u: 

dwc(chi,u) = el( (ke /2 - (e, b - a)) 2 - 2e (e, ka) + kf ). (35) 

It is ease to see that point (0, to) is a singular point on the curve {ch\{u,t)} if 
and only if eo(io) = 0. Therefore the point (0, to) has multiplicity at least two 
as a point of the intersection of two curves {ch\} fl {ch 3 }. 

We will prove that the multiplicity of the intersection of two curves {chi} fl 
{c/13} at the point (0, to) equals to two if eo(to) = 0. For simplicity we assume 
that to = 0. The first equation (|3"3"]) in the local ring R = M.[u, t]( u .t) is 

chi = ci2iu 2 t + CL20U 2 + a\2ut 2 — a\\ut + ao2t 2 , where 

[021, 020, ai2,oii, a a] = [eok(e, a), -fk/2, -el fe/2, e (e, a - 6), 1/2 e§]36) 

and teo = &q. The second equation (|34[) in the local ring i? is 

c/i 2 = u(ch 3 ), ch 3 = &i 2 ut 2 + bunt + b w u + b 02 t 2 + b it, 
where [612, 611, 610, &02, 601] = 

[eofc(e' ; a), -/'e - e e fc(e, a), kf (e' + e' ) , -elk(e', a - b),e e' Q (e, a-b)]. 
An easy computation with MAPLE shows that 

Res t (eft 1, c/i 3 ) = u 2 (K 4 u i + K 3 u 3 + K 2 u 2 + Km + K ). 
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Therefore the point (0,0) has multiplicity two if and only if Kq 7^ 0, i.e. Kq is 
a unit in the local ring R. A straightforward computation shows that 

K = a Q2 (6i 2 a 02 + b 01 2 a 2a + &oi&ioan) = fke^eQ 2 (k(e, e) + ~ (e, a - b) 2 ^j . 

Since /(0) ^ and e£,(0) ^ by the condition we conclude that K ^ 
for a general conic. Hence, the multiplicity i(q, {Res t {hi , h 2 )} fl Q) is equal to 
deg(/ii,t) + 2deg(eo,t) = 4ra. Therefore, we have 

d w (Res t (hi,h 2 )) = 2 deg(Res t {hi, h 2 )) - i(q, {Res t (hi, h 2 )} n Q) 
= 2(5n - 2) - 4n = 6n - 4. 

□ 

Remark 26. We conjecture that the degree of the hypersurfacc is deg V(£)4 
deg(w) — deg(gcd(w)), where w is defined by the formula (|2^|) . 

7 Examples and special cases. 

Let n be the degree of the spine curve £ = {e(t) 6 R 4 }. And let c n {£) be the 
degree of the hypersurface T{£ ) with the spine curve £. 

Polynomial case. Assume that the spine curve is polynomial, i.e. eo = 1. 
By the theorem in Xu et al.(2006)| the degree of hypersurface T(£) with the 



polynomial spine is at most 4n — 2. 

For the general spine curve we have c n {£) = 6n — 4, i.e. c n (£) < dm — 4. The 
lower bound is not clear. There are examples of spine curves with the following 
degrees: 

c 2 {£) = 3,4,5,6,8; 

c 3 (£) = 6,7,8,9,10,11,12,14; 

c 4 (£) = 8,9,10,11,12,13,14,15,16,17,18,20. 

It seems that there does not exist a spine curve such that c n (£ ) = 6n — 5. 
We consider three examples. 



Example 27. Let us consider the following spine curve: e(t) = (0, 0, j^s,^: 



3// 



l+t 

n = 2. This is a proper parametrization of an ellipse in R . We find the Plucker 
coordinate vector P = £{t)C A £'{t)C and q — gcd(P) = 1. Also, we see that 
degV(£) = deg P — deg(q,t) = 8 and 7 = deg(w) = 2. If we run the /x-basis 
algorithm with two input vectors £(t)C,£'(t)C we get the output two vectors 
Ei and E 2 : 

E 1 -y T = 4t 3 y 3 + {-u-41y )t 2 + 12ty 3 + 9y Q -u-6y 4 , 
E 2 -f = (w-9yo-62/4)i 3 -m 2 2/3 + (4l2/o + w)t-4y 3 - 

Now we can find the implicit equation G = Res (Si • y T ,E 2 ■ y T ,t) of the dual 
variety V{£ ). The polynomial G contains 26 monomials and has degree 6 (as 
in Proposition [53]) . The equation of the hypersurface T{£) is defined by the 
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polynomial 

F{yo,-,V4.) = y 2 G{(y,y)/y Q ,y , Vi, 2/2,2/3,2/4) of degree 8. Since, 

F(l, 2/1,2/2, J/3,0) = (2/i 2 + 16 + 2/2 2 + 8 2/ 3 + 2/3 2 ) (2/1 2 + 2/2 2 + 16 - 8 y 3 + y 3 2 ) 

(-225 + 25 2/i 2 + 25 2/ 2 2 + 9 2/3 2 ) 2 , 

the 0-envelope of the canal surface Envo(£) = T(£ ) n {2/4 = 0} is reducible. The 
canal surface C is the double ellipsoid of revolution (—225 + 25?/i 2 + 25?/2 2 + 
9 J/3 2 ) 2 ■ Indeed, for the computation of C(£) we should assume that the variable 
2/4 = and to repeat the same steps as above. We should consider only the 
first 5 coordinates of the vectors £(t)C,£'(t)C. Let us denote these two vectors 
with 5 coordinates by Z?i,£>2- But this time we see that the Phicker vector 
P = D\ A D 2 has a non-trivial common divisor, i.e. q = gcd(P) = t 2 — 1. So, 
using the /x-basis algorithm we find the /x-basis i?i,i?2 for the input D\,D 2 - 
In this case we see that degi?i = degi?2 = 2. Now we find the resultant G = 

Res t (Ri-y T ,R2-y T ) = (l6j/ 3 2 + 225?/o 2 - 25j/ -u) , wherey = (u, 2/0,2/1,2/2,2/3) ■ 
After the substitution u = (yf + y 2 + y 2 )/yo we obtain the implicit equation 
of the canal surface the double ellipsoid (—225 -I- 25 2/1 2 + 25?/2 2 + 9?/3 2 ) 2 - We 
can see this geometrically, too. The point e(t) G K 4 corresponds to the sphere 
S(e(t)) G K 3 with a center on the 2/3-axis. If t G [—1/2,1/2] then the sphere 
S(e(t)) is tangent to the ellipsoid EL = (—225+25 yi 2 +25 1/ 2 2 +9 y 3 2 ), and inside 
this ellipsoid. Moreover, the real envelope of the family S(e(t)),t G [—1/2, 1/2] 
is the ellipsoid EL. Note that the sphere S(e(l/t)) has the same center but 
the opposite radius to the sphere S(e(t)), i.e. it has the opposite orientation. 
Therefore, the real envelope of the family S(e(t)), t G (—00, —2] n [2, 00) is the 
same ellipsoid EL. Hence, from the point of Laguerre geometry the envelope of 
the whole family S(e(t)) is the double ellipsoid EL 2 . Note, that the d-offset to 
the canal surface, in this case is the d-offset to ellipsoid EL and it has degree 8. 
Also, we can check that by Theorem the degree of the T(£) hypersurface is 8, 
too. For a detailed study and other examples of canal surfaces with a quadratic 
spine curve we recommend to look at the paper [Krasauskas, Zube(2007)| . 

Example 28. Consider the polynomial spine curve e(t) — {it 2 + 1, At 2 + t, 0, 5t 2 ) , 
n = 2. We find the Phicker coordinate vector P = £(t)C A £'(t)C and q = 
gcd(P) = 1. Also, we see that deg V(£ ) — degP—degq — 4 and7 = deg(w) = 1. 
If we run the /i-basis algorithm with two input vectors £(t)C, £'{t)C we get the 
output of two vectors 

E x = -u+ (-l-7t 2 -8t 3 )y + (2 + 6t 2 ) yi + 2t{l + 4t)y 2 ~l0t 2 y 4 , 
E-2 = -<(7+12<)y + 6%i + (l + 8t)2/2-10<2/4, 

and find the implicit equation G of the dual variety V(£) (it contains 54 mono- 
mials, so we do not present an explicit formula). Finally, we find that 02(8) = 
d w {G) = 5. For this example, we have degC = degT(£), i.e. the implicit 
degree of the canal surface is 5. Note that this contradicts Theorem 4 in 
[Xu et al. (2006)| , because in this example the degree of the canal surface is an 
odd number. It seems that the mentioned theorem gives only an upper bound 
estimation, but not the exact degree formula of canal surfaces with polynomial 
spine curve. 
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Example 29. In the next example we take the following spine curve: 

/ (i-t 2 ) 2 t(i-t 2 ) \ 
e M = ( n+t 2 ) 2 ' ^ (l+Pr ' ^ TT* 7 ' ^ ) ' n = 4- The first three coordinates define 

the Viviani curve, i.e. it is intersection curve of the sphere and the tangent 
cylinder. We find the Plucker coordinate vector P = £(t)C A £'(t)C and q = 
gcd(P) = 1. Also, we see that degV(£) = degP — degg = 6. If we run the 
/^-basis algorithm with two input vectors £(t)C,£'(t)C we get output of two 
vectors 

E x = (0,4 + 4i 2 ,4-4t 2 ,6t-2i 3 ,6t + 2t 3 ,4 + 4i 2 ) , 

E 2 = (0,4i + 4i 3 ,4i(-l + t 2 ) ,2-6i 2 ,2 + 6£ 2 ,4£ + 4t 3 ) , 

both of degree 3 and find the implicit equation G of the dual variety V(£) (it 
contains 58 monomials). Finally, we find that Ci(£) = d w (G) = 10. For this 
example, we have degC = degT(£), i.e. the implicit degree of the canal surface 
is 10. 
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